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We analyse the CP asymmetry of the B — > (f>K s and B — > rfKs processes in general supersymmetric models. We show that chromomag- 
netic type of operator may play an important role in accounting for the deviation of the mixing CP asymmetry between B — > (j)K s and 
B — > J/if/Ks processes observed by Belle and BABAR experiments. We also show that due to the different parity in the final states of 
these processes, their supersymmetric contributions from the R-sector have an opposite sign, which naturally explain the large deviation 
between their asymmetries. 
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1 Introduction 

One of the most important tasks for B factory experiments 
would be to test the Kobayashi-Maskawa (KM) ansatz for 
the flavor CP violation. The flavor CP violation has been 
studied quite a while, however, it is still one of the least 
tested aspect in the standard model (SM). Although it is 
unlikely that the SM provides the complete description of 
CP violation in nature (e.g. Baryon asymmetry in the uni- 
verse), it is also very difficult to include any additional 
sources of CP violation beyond the phase in the CKM mix- 
ing matrix. Stringent constraints on these phases are usu- 
ally obtained from the experimental bounds on the electric 
dipole moment (EDM) of the neutron, electron and mer- 
cury atom. Therefore, it remains a challenge for any new 
physics beyond the SM to give a new source of CP vio- 
lation that may explain possible deviations from the SM 
results and also avoid overproduction of the EDMs. In su- 
persymmetric theories, it has been emphasised [ 1 1 that 
there are attractive scenarios where the EDM problem is 
solved and genuine SUSY CP violating effects are found. 

Recently, BABAR and Belle collaborations announced a 
2.1 a deviation from sin 2/3 in the B — » cpK s process [|3|3|. 
In the SM, the decay process of B — > (pK is dominated by 
the top quark intermediated penguin diagram, which do not 
include any CP violating phase. Therefore, the CP asym- 
metry of B — > J/ifsKs and B — > cpKs in SM are caused 
only by the phase in B° - B mixing diagram and we ex- 
pect S j/^Ks = S<pK s where S f CP represents the mixing CP 
asymmetry. The B — > rj'Ks process is induced by more dia- 
grams since r{ meson contains not only ss state but also uu 
and d d states with the pseudoscalar mixing angle P . Nev- 
ertheless, under an assumption that its tree diagram con- 
tribution is very small, which is indeed the case, one can 
expect = Srj'Kg [ 2 4| as well. Thus, the series of new 
experimental data surprised us: 
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It was pointed out [ 5 ] that the discrepancy between Eq. |T} 
and Eq. might be explained by new physics contribu- 
tion through the penguin diagram to B — > (pKg . However, 
in that case, a simultaneous explanation for the discrepancy 
between Eq. (0 and Eq. (0 is also necessary. We show 
our attempts to understand all the above experimental data 
within the Supersymmetric models. 

2 The mass insertion approximation 

As mentioned, the SUSY extension of the SM may provide 
considerable effects to the CP violation observables since 
it contains new CP violating phases and also new flavour 
structures. Thus, SUSY is a natural candidate to resolve the 
discrepancy among the observed mixing CP asymmetries 
in 5-meson decays. 

In the following, we will perform a model independent 
analysis by using the mass insertion approximation [|6|. 
We start with the minimal supersymmetric standard model 
(MSSM), where a minimal number of super- fields is intro- 
duced and R parity is conserved, with the following soft 
SUSY breaking terms 
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- ^(m 3 §g + m 2 W a W a + miBB^ , (4) 



where i, j are family indices, a,b are SU (2) indices, and 
€ab is the 2 x 2 fully antisymmetric tensor, with en = 1. 
Moreover, (p a denotes all the scalar fields of the theory. Al- 
though in general the parameters yu, B, A a and m; can be 
complex, two of their phases can be rotated away. 
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The mass insertion approximation is a technique which is 
developed to include the soft SUSY breaking term with- 
out specifying the models in behind. In this approxima- 
tion, one adopts a basis where the couplings of the fermion 
and sfermion are flavour diagonal, leaving all the sources 
of flavour violation inside the off-diagonal terms of the 
sfermion mass matrix. These terms are denoted by (A^ g ) u , 
where A, B = (L,R) and q = u,d. The sfermion propagator 
is then expanded as 



1 Sab 



k 2 - m 2 (k 2 - in 2 ) 2 



(5) 



where 1 is the unit matrix and m is the average squark 
mass. The SUSY contributions are parameterised in terms 
of the dimensionless parameters (S q AB )ij = (A AB ) lj /m 2 . This 
method allows to parametrise, in a model independent way, 
the main sources of flavor violations in SUSY models. 

Including the SUSY contribution, the effective Hamilto- 
nian for the penguin diagrams are written as 
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where C 3 ~ C 6 (C 3 ~ C 6 ) include (S LL ) 23 ((S RR ) 23 ) contribu- 
tions and C g {C g ) include (S ll , lr ) 23 ((Srr, rl ) 23 ). The terms 



with tilde are obtained from Q 
L^R. 



and Oi g by exchanging 



3 Can we explain the experimental data of 

S^ Ks in SUSY? 

Following the parametrisation of the SM and SUSY ampli- 
tudes in Ref.[0, S ^ can be written as 



>4>Ks 



sin 2/3 + 2R<f, cos S 12 sin(# + 2/3) + R% sin(2# + 2/3) 
1 + 2R(p cos 5\ 2 cos 0(p + R^ 

where R<p = |A SUSY /A SM |, = arg(A sus Y / A SM ) , and S 12 is the 
strong phase. 

We will discuss in the following whether the SUSY contri- 
butions can make S<f,K s negative. Note that the mass inser- 
tions {Sab)i3 have already been constrained by the experi- 
mental data for Br(B — > X s y) : 



\(S L L,Rr)23\ < 1, |(<W.)23l < 1.6 X 10 



For Mc, - m~ g = 500 GeV, we obtain 
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A SM 



: 0.23(^)23+0.23(^)23+97.4(^)23+97.4(^)23(10) 



The constrains from Br(B — > X s y) gives the maximum 
|A SUSY |/|A SM |: 
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Figure 1. (a)0 3 - 6 contributions which include (S LL , RR )23 mass 
insertions, (b) O g contribution which includes (Sll,rr,lr,Rl)23 
mass insertions. 



As emphasised in [ 7 1, the leading contribution to AB = 1 
processes come from the chromomagnetic penguin opera- 
tor Og(O g ), in particular from the part proportional to the 
LR (RL) mass insertions which is enhanced by a factor 
rrig/mt,, where C g C g are given by 
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Note that the mass insertions appearing in the box diagrams 
are (Sab)\3 (A,B = L or R), thus, SUSY contributions 
to box diagram and to penguin diagram are independent. 
S ^ sin 2/3 indicates the smallness of (Sab)u [ 8|. 
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< 1.6 For (S LR , RL ) 23 
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In Fig 13 we present plots for the phase of (S d LL ^ RK) ) 23 and 
(S d LR ^ RL) ) 23 versus the mixing CP asymmetry S^k s when the 
strong phases are ignored. We choose the three values of 
the magnitude of these mass insertions within the bounds 
from the experimental limits from B — > X s y. Each plot 
shows a contribution from an individual mass insertion by 
setting the other three to be zero. As can be seen from these 
plots, the LR (same for RL) gives the largest contribution 
to S (pK s • I n order to have a sizable effect from the LL or 
RR, the magnitude of (S d LL{RK) ) 23 has to be of order one and 
furthermore, the imaginary part needs to be as large as the 
real part. In any case, it is very difficult to give negative 
value of S(f,K s from (S d L )i3 or (S RR ) 23 mass insertion. If 
the experimental data remains as small as the current val- 
ues, LL, RR dominated models would get sever constrains 
on some parameters. Note that S^k s decreases as SUSY 
masses becomes smaller. In [ 9|, a choice of niq ^ 350 
GeV has been used and a negative S ^ for LL, RR models 
has been obtained. 
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Since the coefficient for each mass insertions are similar, 
we use the following definition to simplify our following 
discussions: 



iSUSY 



= St + S R 



4>K S 



tSUSY 



A SM 



(16) 



rfK s 



Figure 2. Result for S ^ Ks in terms of the phase in mass insertions. 

4 What happened to the B — > rfK s process? 

Although B —> (pKs and B —> rfKs are very similar pro- 
cesses, the parity of the final states can deviate the result. In 
the naive factorisation approximation, the amplitudes are 
written as a product of Wilson coefficients, form factors 
and decay constants: 



A(B -> (f>(?]')K) oc C Wilson F 



(13) 



The decay constants appear in the calculation by sandwich- 
ing the V+A current (Of and 0[ contributions, respectively) 
with (pirf) and vacuum: 



(0|^(1 ± y 5 )s\(f>) = m^ffy 
(01^(1 + 75 )s\rj') = ±ifr 1 >P» 



(14) 
(15) 



As can be seen from Eqs. JlH and fT5l . the vector meson 
(p picks up the y^ term while the pseudoscalar meson r{ 
picks up the y^ys term so that contributions from 0[ and 
0[ obtain opposite signs for . 



where S L includes contributions from (S LL )23 and (6 lr) 23 
and 6r includes contributions from (S RR )23 and {8rl)i3- 
Now let us show how this sign flip effects to the mixing CP 
violation S^ks and S^Ks- Since both Sl and S R are com- 
plex number, we have four parameters to be fixed while we 
have only two experimental data. Thus, we fix two param- 
eters and perform a case-by-case study in the following. 



• Case 1: \6 R \ » \5 L 
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• Case 2: \S L \ = \S R \ (AO = arg£ L - zrg6 R ) 
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• Case 3: arg5 L = argfc (A\S\ = \S L \ - \8 R \) 
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Figure 3. Schematically described naive factorisation approxi- 
mation for B — > (pK s and B — > tj'K s processes. 



As a result, the sign of the RR and RL contributions are 
different for B <f>K s and B -» rfK s f fTOl : 



0.23(^)23+0.23(^)23+97.4(^)23+97.4(^)23 



<PK S 



[-^T) = 0.23(^)23-0.23(^)23 + 101(^)23-101(^)23 
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QS L \ + \6r\W"* 6l 
A\S\e iarg6L 



• Case 4: argfo = argS L + n/2 (tana = \6 R \/\6 L \) 



iSUSY 



<PK S 
^SUSY \ 



Vfcl 2 + \6 R \ 2 e^ +a) 
Vfcl 2 + \S R \ 2 e^- a) 



In Fig. 01 we show some examples of the parameter sets 
with which we can reproduce both experimental data of 

S <f,K s and SrfKs- 
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0<pK s 



Figure 4. Case 1: \6 R \ dominating. Case 2: \6 R \ = \6 L \ with 
A0 = tt/10. Case 3: avg6 L = arg6 R with A\6\ = +0.2. Case 
4: arg£# = arg£ L + tt/2 with or = 3tt/8. Solid line: 5^ for 
|A SUSY |/|A SM | = 0.5. 

5 On the branching ratio of B ^> rj'K s : Glu- 
onium vs. New physics 

In 1997, CLEO collaboration reported an unexpectedly 
large branching ratio lITTl 

Br^ (B° -> K°rf) = (89!}* + 9) X 10 -6 (17) 

which is confirmed by Belle [[12| and BABAR [0): 

BELLE = (19 + _\ 2 6 + 8) x 10" 6 , (18) 
BABAR = (76.9 + 3.5 + 4.4) x 10 -6 (19) 

Considering the theoretical prediction by the naive factori- 
sation approximation 

Br theo (B -> Krf) ^ 25 x 10 -6 , (20) 

the experimental data is about factor of three large, thus, 
there have been various efforts to explain this puzzle. On 
one hand, new physics contributions have been discussed 
[ 13 1 . However, the enhancement by new physics con- 
tributions through penguin diagrams ends up with large 
branching ratios for all other penguin dominated processes. 
Therefore, one needs a careful treatment to enhance only 
B — > v[K process without changing the predictions for the 
other processes. On the other hand, since this kind of large 
branching ratio is observed only in B — > r(K process, the 
gluonium contributors which only exist in this process have 
been a very interesting candidate to solve the puzzle [ [14| 
[ 15 1 though the amount of gluonium in v[ is not precisely 
known [ 16 1. In this section, let us discuss the effect of our 
including SUSY contributions to the branching ratios for 
B^(f>K<mdB^ r/K. 

Inclusion of the SUSY contributions modify the branching 
ratio as: 

£ r s M + susY = Br su x fl + 2 cos Si]SY R + R 2 ] 



where R - |A SUSY |/|A SM |. As we have shown, to achieve 
a negative value of S^, we need # SUSY ^ -tt/2, which 
suppresses the leading SUSY contribution. As a result, for 
instance, R = 0.5 leads to: 

Br(B -> (f>K s ) = (7.8 x 10" 6 ) x 1.25 = 9.7 x 10" 6 (21) 

which is within the experimental data (9. 1 + 2.6) x 10" 6 . On 
the other hand, the phase for B — > rj'K is different from the 
one for cpKs , as is discussed in the previous section. For 
instance, Case 2 gives us the maximum value of: 

Br(B -> rj'Ks) = (25 x 10 -6 ) x 3.25 = 81 x 10 -6 (22) 

However, this kind of enhancement would appear all the 
other two pseudo-scalars channels (such as B — > Kn) and 
might cause some problems. 

As a whole, we would like to suggest that the solution for 
the branching ratio puzzle is not only the SUSY contribu- 
tion but a combination of SUSY contribution and the glu- 
onium contribution. Here, let us show the dependence of 
the gluonium contribution to the S^k s - Including the glu- 
onium contribution (see Fig. |5j, the amplitude is modified 
to 

A = A™ Ks + A*Y£ + G SM + G SUSY (23) 

where G SM and G SUSY are the new mechanism contributions 
to SM and SUSY, respectively. Let us parametrise the un- 
known gluonium content in v[ as r = G SM /A S ™ K . Our result 
is shown in Fig. [6|when we vary r from to 0.3. As can be 
seen from this figure, the dependence of S^k s on r is not 
very strong, therefore, we can enhance the branching ratio 
by gluonium contribution without disturbing our findings 
for SrfK s m me previous section. 




Figure 5. A contribution from gluonium content in 77' to the B — > 

rj'K process. 
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Figure 6. Branching ratio versus the mixing CP asymmetry in the 
B — > 77^ process. The parameter r represents the contribution 
from the gluonium diagram. 

6 Conclusions 

We studied the supersymmetric contributions to the CP 
asymmetry of B — > (pKg and B — > tj'Ks in a model inde- 
pendent way. We found that the observed large discrepancy 
between S j/^k s an d can be explained within some 
SUSY models with large (#l/?)23 or (Srl)23 mass insertions. 
We showed that the SUSY contributions of (Srr)23 and 
(^rl)23 to B — > 0^ and 5 — > 77' £5 have different signs. 
Therefore, the current observation, < S^k s , favours 
the (Srr,rl)23 dominated models. We also discussed the 
SUSY contributions to the branching ratios. We showed 
that negative S^ Ks and small SUSY effect to Br(B -> cpK) 
can be simultaneously achieved. On the other hand, we 
showed that SUSY contribution itself does not solve the 
puzzle of the large branching ratio of B —> tj'Ks. We in- 
cluded the gluonium contributions to B — > tj'Ks • We found 
that our conclusion for S ^k s does not disturbed by gluo- 
nium contributions. As soon as the experimental errors are 
reduced, the CP violation of B — > cpKs and B — > tj'Ks 
will be able to give a strong constraints on the (Sab)23 mass 
insertions. 
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